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We investigate the direct production of supersymmetric scalar taus at the LHC. We present the 
general calculation of the dominant cross section contributions for hadronic stau pair produc- 
tion within the MSSM, taking into account left-right mixing of the stau eigenstates. We find 
that /7-quark annihilation and gluon fusion can enhance the cross sections by more than one or- 
der of magnitude with respect to the Drell-Yan predictions. For long-lived staus, we consider 
CMSSM parameter regions with such enhanced cross sections and possible consequences from 
recent searches. We find that regions of exceptionally small stau yields, favoured by cosmology, 
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1. Introduction 

The ongoing experiments at the Large Hadron Collider (LHC) set limits for new physics at the 
TeV scale in an unprecedented way. Amongst the most promising candidates for new physics are 
supersymmetric (SUSY) theories. Within this class of theories, scenarios with a scalar tau z being 
the lightest ordinary SUSY particle (LOSP) are often considered. In R-parity-conserving scenarios 
the lightest SUSY particle (LSP) is then usually assumed to be extremely weakly interacting, e.g., 
a gravitino or an axino. Such scenarios offer spectacular phenomenological signatures at colliders: 
long-lived charged massive particles (CHAMPs). Furthermore, they can be associated with in- 
triguing early Universe cosmology, cf. [2] and references therein. For example, an overly abundant 
charged relic can spoil the successful predictions of BBN. Resulting strong bounds can be evaded 
in parameter regions where the stau annihilates efficiently, e.g., via a heavy Higgs resonance or due 
to large stau-higgs couplings [3,4]. In scenarios with broken R-parity, the f can also be the LSP 
that decays promptly (or delayed in scenarios with partly broken R-parity) into Standard Model 
(SM) particles. 

In this work we present theoretical predictions for direct f j f j* production at the LHC within the 
minimal SUSY extension of the Standard Model (MSSM) independent of a possible longevity of 
the stau. We include Drell-Yan processes as well as ^-quark annhilation and gluon fusion. Special 
attention is given to mixing effects due to large Yukawa couplings. After evaluating numerical 
results in the phenomenological MSSM to investigate possible enhancement effects, we interpret 
our results in the CMSSM assuming a long-lived fi. This proceeding is mainly based on Ref. [1] 
where also more references can be found. One new aspect which has not been discussed in [1] is 
the interpretation of recent experimental exclusion limits set by the CMS experiment at the LHC. 

2. Direct production of stau pairs at the LHC 

Within the MSSM, stau pairs can be produced directly at the LHC, 

PP^iji*, (2.1) 

where f, ; denotes any of the two stau mass eigenstates. Due to off-diagonal elements in the stau 
mass matrix large mixing between gauge eigenstates (Tlj^r) can occur for the mass eigenstates 
(fi , T2) and may not be neglected for third generation sleptons. This mixing can become large and 
is proportional to the SUSY parameters pL, tan j3 and A x . Here we concentrate on the direct produc- 
tion of the lighter f 1 fj* pairs, particularly in parameter regions with large mixing. Production rates 
of f 2 f| and similarly f l f| can be obtained in close analogy. 

At the LHC the leading contribution to direct i l fj* production up to order &{a 2 a 2 ) are given 
by the following channels: 

• qq induced Drell-Yan type processes at @(a 2 \ see figure |l](a), and corresponding NLO 
(SUSY)QCD corrections at @(a 2 s a 2 ). 

• bb annihilation, mediated by the neutral gauge bosons (7, Z) and by the neutral CP-even 
Higgs bosons (h°, H°) at &{a 2 ) (however suppressed by the low bottom-quark PDFs), shown 
in figure [j](b). 
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Figure 1: Feynman diagrams for stau pair production (a) via the Drell-Yan process and (b) via bb annihila- 
tion. Here q = u,d,c,s. 

9 9 

9 9 9 * 

Figure 2: Feynman diagrams for the gluon fusion contribution to stau pair production. The quarks q and 
squarks ; = 1,2, running in the loops can be of any flavor. 

• gluon-gluon fusion at order 0(aja 2 ), mediated by a quark or squark loop, as shown in 
figure ^. 

The Drell-Yan production cross section at leading order depends only on the stau mass 
and the stau mixing angle 6? . The parametric dependence of the NLO corrections on the SUSY 
parameters in general is small. Although suppressed by the low bottom-quark density inside pro- 
tons the bb channel can be enhanced by on-shell Higgs propagators and by the bottom-Higgs and 
the stau-Higgs couplings in certain regions of the SUSY parameter space. In the basis of mass 
eigenstates, the stau-Higgs coupling is proportional to sin 6? and thus becomes important for large 
mixing in the stau sector. Both the stau-Higgs and the bottom-Higgs couplings are proportional 
to jU and the respective trilinear coupling A % / b . In order to study possible enhancement effects in 
direct f j production due to large stau-Higgs couplings, parameter regions with relatively large 
tan/3 are to be considered. Here radiative corrections to the bbh /bbH° vertex can be important and 
drive down the cross section compared to the tree-level result. Leading tan j3 -enhanced corrections 
can be resummed to all orders in pertubation theory by using an appropriate effective bottom-quark 
mass, m| ff , and effective bbh Q /bbH° couplings. Here we adopt this approach, as explained in detail 
in appendix B of [1]. Finally, even though gluon-induced contributions are formally of higher or- 
ders, tff(aja 2 ), they can give sizeable contributions at the /jp-machine LHC at high center-of-mass 
energies where the gg luminosity is significantly higher than the qq luminosity. As the bb channel, 
the gluon-gluon channel depends strongly on the parameters in the stau-Higgs sector. 

We use the programs FeynArts 3 . 6 and FormCalc 7 . with LoopTools 2 . 6 to 
generate and calculate the amplitudes corresponding to the Feynman diagrams of figures [I] and |[ 
The Higgs boson masses and the H° width are computed with FeynHiggs 2.7.4. The QCD 
and SUSY-QCD corrections at NLO for the Drell-Yan channel calculated in [5] we evaluate with 
Prospino 2 and scale our cross sections with the resulting ^-factors. Furthermore, we use a re- 
summed effective bbh° /bbH° vertex for the gluon fusion and bb contributions, as explained above. 
Since we do not include (expected to be positive) higher-order QCD and SUSY-QCD corrections to 
the Higgs-mediated channels, our analysis gives a conservative estimate of the enhancement effects 
from the to-annhilation and gluon fusion production channels. 
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At the LHC, the gluon-fusion and &Z?-annhilation processes with an s-channel Higgs boson can 
become resonant in regions of the SUSY parameter space in which the Higgs boson H is heavier 
than the two produced staus. In parameter regions with m H o > 2m^ we therefore include the total 
decay width of the H° boson, in the propagator, 

1 1 

P 2 ~ m 2 H o P 2 ~ m 2 H o + im H oT H o ' 

3. Numerical results 

The cross section for direct stau production depends mainly on m^, m H o, tan/3, and on df (or 
equivalently on /J. and A T ). Let us now investigate the dependence on these parameters. 

As a starting point, we choose a fpLOSP scenario with moderate squark masses and a large 
stau-Higgs coupling, fixed by the following soft-breaking parameters at the low scale: 

M l =M 2 =M 3 = 1.2 TeV, A t = A b = A x = 600 GeV, 

(3.1) 

m Qi = m U, = m D, = 1 TeV > m L 1/2 = m E l/2 = 500 GeV > 

If not otherwise stated, we choose, 

i= 45°, ra f ,=200GeV, 

(3.2) 

tan p = 30, pi = 500 GeV, m A = 400 GeV, 

as inputs for the third-generation sleptons and the Higgs sector. 

From these input parameters, we calculate the physical MSSM parameters using tree-level re- 
lations for sfermions, neutralinos, and charginos. Physical masses are then passed to P ro spino 2 
to calculate the Drell-Yan K-factors at NLO in QCD and SUSY-QCD. The NLO corrections to the 
Drell-Yan channel typically amount to 20 — 40 % in the considered parameter space. 

In figure || we show the direct production cross section for f i -pairs at the LHC with \/S = 
14 TeV as a function of (a) m^, (b) di , (c) m H o, and (d) tan/3. In figures |3|(b) and (d), we move to 
mf, = 190 GeV where stau production is possible via an on-shell H°. The dashed (red) lines show 
the Drell-Yan (DY) cross section at NLO, whereas the solid (blue) lines include the additional bb 
and gg contributions. The Drell-Yan cross section depends on and Of only. It decreases strongly 
for increasing fi masses and varies roughly by a factor slightly larger than 2 with df , as shown in 
figure |3|(b), being largest for 0, i.e., an almost left-handed ?i. 

The impact of the bb and gg channels depends strongly on the mass hierarchy between f i and 
H°, as can clearly be seen in Figures ||(a) and (c). If m H o > 2m^, these additional channels can 
change the direct production cross section by more than one order of magnitude with respect to the 
Drell-Yan result. At the threshold m H o = 2/nf, , the bb and gg contributions drop steeply and are 
only marginally important for m H o <C 2mf, . 

Figures ||(b) and (d) illustrate the dependence of the total direct production cross section on the 
parameters and tan /3 that govern the stau-Higgs-coupling strength. The additional contributions 
from the bb and gg channels are tiny in cases of very small mixing, 6? — > 0, %, but become most 
important for maximal mixing, i.e. , at %/A. For very large tan /3, additional contributions push 
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Figure 3: Cross section of direct fi production at the LHC with y/S = 14 TeV as a function of (a) m{,, 



(b) Of, (c) m H o, and (d) tanj3. No kinematical cuts are applied. SUSY input parameters as given in (3.1) and 



(3.2). Note that here m H o rs Ma (decoupling limit). 



up the total direct production cross section by up to two orders of magnitude and are still sizeable 
for small tan/3. 

Let us now turn to a scenario where the H° is very heavy and thus almost decoupled, m H o = 
1 TeV. We again investigate the dependence of the total cross section on df and tan /3 , shown in 
figure |], where we focus on enhanced fi-fj"-^ couplings. In figure f|(a) we choose jj, = 800 GeV 
and tan/3 = 50 to increase this coupling. All other parameters are fixed according to (3J_) and 



(p.2|). Again, contributions from additional bb and gg channels can be sizeable and the enhancement 
amounts to a factor of two or three when considering very large values of tan /3 and maximal mixing 
6f « 71 / 4. However, here dominant contributions come mainly from off-shell h diagrams together 
with large couplings. Thus, the relative importance between the gg channel and the bb channel 
can be different compared to on-shell H° production, as the two Higgses couple differently to the 
squark loops. We want to note that, despite the large couplings, all considered parameter points 
are in agreement with bounds from the potential occurrence of (color and) charge breaking (CCB) 
minima [6]. 
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Figure 4: Cross section of direct X\ production at the LHC with \/S = 14 TeV as a function of (a) Of and 



(b) tan/3. No kinematical cuts are applied. SUSY input parameters as given in (3.1) and (3.2) if not stated 
otherwise. Note that here m H o f=s = 1 TeV (decoupling limit). 




We summarize the potential impact of the bb and gg channels again in figure ||, where the 
Drell-Yan contribution (dashed lines) and the full cross sections (solid lines) are shown for f j fj* 
production at the LHC for y/S = 14 TeV (top, red) and for y/S = 1 TeV (bottom, blue). When 
going down from 14 TeV to 7 TeV, the cross section decreases by up to about a factor of 5. The 
relative contribution of the bb and gg channels however can be similarly important. Thus, for both 
y/S = 7 TeV and 14 TeV (and also = 8 TeV), the bb and gg channels should not be neglected in 
a precise cross section prediction. 
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Figure 6: Colored contours in the CMSSM wo-mi/2 plane with tan/3 = 55, Aq = 2mo, and fx > of the 
total direct T] f* production cross section (in femtobarn) after cuts at the LHC with y/S = 7 TeV. The white 
region in the lower left is excluded because of a tachyonic spectrum, impossible EWSB, or < 82 GeV. 
In the upper white area the lightest neutralino %\ is the LOSP. The lower hatched area is in tension with 
CCB constraints. The solid black lines indicate the mass of the lighter Higgs boson h°. The dashed black 
lines are contours of constant X\ masses, the white lines refer to the heavier Higgs boson H°. The hatched 
area at small f i masses is excluded by recent CMS searches. 

4. Direct production of long-lived staus within the CMSSM 

Now we turn to the CMSSM as a benchmark model. We consider the mo-mm plane of the 
CMSSM with Aq = 2m , tan/3 = 55, and }X > 0. Here one often finds a t\ LOSP. Moreover, 
this parameter choice is cosmologically motivated by the possibility of a long-lived fi, due to an 
extremely weakly interacting dark matter candidate such as the axino or gravitino and by excep- 
tionally small stau yields [4]. In this plane, tan/3 is large and the Z\ prefers to be right-handed 
(8f > 71 / 4; this is generic in the CMSSM due to the different running of the left-handed and right- 
handed soft masses), and the bb and gg channels can give large contributions to the stau production 
cross section. In the following, our cross section predictions include the Drell-Yan channels with 
NLO K-factors, the bb annihilation and the gg fusion contributions. Additionally we apply the 
following experimental cuts used by collider experiments to discriminate from SM backgrounds in 
their CHAMP searches: 

/? T >40GeV, 0.4</3<0.9, |t]|<2.4, (4.1) 

where p T denotes the transverse momentum, /3 the velocity, and v\ the pseudorapidity of the pro- 
duced CHAMP. These cuts reduce the cross section by about 30%. 
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The resulting cross section is shown in figure ^ It depends mainly on and m H o and varies 
over several orders of magnitude in the given parameter range. At the LHC with y/S = 7 TeV, it 
reaches 10 2 fb for mo < 700 GeV, < 500 GeV and drops to 10~ 3 fb for, e.g., mo ~ m X j 2 ~ 
2 TeV. 

In a recent analyses [7] using an integrated luminosity of Jzf = 4.7fb~ 1 , the CMS experiment 
excludes long-lived staus with 

m fl < 232GeV. (4.2) 

A thereby excluded parameter region of the considered CMSSM plane is indicated in figure ||. 
Comparing with figure 12 of [1], this new limit seems to disfavour the cosmologically appealing 
possibility of exceptionally small stau yields in this CMSSM plane (and also in the CMSSM in 
general). The limit ( |4~2| ) set by the CMS collaboration is based on the assumption of directly pro- 
duced %\ pairs via the Drell-Yan channel at NLO and should thus be considered as fairly model 
independent. Determination of the exact discovery reach and/or the exclusion limits including the 
additional bb and gg channels should be performed in context of a detailed study including detector 
effects. 

5. Conclusions 

We have studied the direct hadronic production of a pair of staus fifj* within the MSSM. In 
addition to the well-known Drell-Yan process, we considered production processes initiated by 
bb annihilation and gluon fusion, with all third-generation mixing effects taken into account. In 
parameter regions with non-negligible mixing these contributions can enhance the direct production 
cross section significantly and should always be included. These predictions are independent of the 
stau lifetime and applicable in %\ LSP scenarios with the fi being the next-to-LSP as well as in 
settings in which the fi is long-lived. 

Within the CMSSM and assuming the fi being long-lived on the scale of colliders, we have 
provided cross section predictions for direct stau production at the LHC with \f~S = 1 TeV. Recent 
exclusion limits from the experiments at the LHC are interpreted in this framework. From this, at 
least within the CMSSM, the possibility of exceptionally small stau yields seems to be disfavoured. 
The current \/S = 8 TeV run of the LHC will fully probe the viability of such yields within the 
CMSSM. 

Once a deviation from SM backgrounds is observed, the parameters of the underlying theory 
have to be determined. In our detailed analysis [1] we offer various ideas how, here, direct stau pro- 
duction can be of utter importance. Measuring the enhancement due to bb annihilation and gluon 
fusion together with kinematic distributions can even offer the possibility to test early Universe 
implications in the laboratory. 
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